Miscibility of Phosphatidylcholines in Bilayers: Effect of Acyl Chain Unsaturation

The miscibility of phospholipids in a hydrated bilayer is an issue of fundamental importance for understanding the organization of biological membranes. Despite research on lipid miscibility, its molecular basis remains poorly understood. In this study, all-atom MD simulations complemented by Langmuir monolayer and DSC experiments have been performed to investigate the molecular organization and properties of lipid bilayers composed of phosphatidylcholines with saturated (palmitoyl, DPPC) and unsaturated (oleoyl, DOPC) acyl chains. The experimental results showed that the DOPC/DPPC bilayers are systems exhibiting a very limited miscibility (strongly positive values of excess free energy of mixing) at temperatures below the DPPC phase transition. The excess free energy of mixing is divided into an entropic component, related to the ordering of the acyl chains, and an enthalpic component, resulting from the mainly electrostatic interactions between the headgroups of lipids. MD simulations showed that the electrostatic interactions for lipid like-pairs are much stronger than that for mixed pairs and temperature has only a slight influence on these interactions. On the contrary, the entropic component increases strongly with increasing temperature, due to the freeing of rotation of acyl chains. Therefore, the miscibility of phospholipids with different saturations of acyl chains is an entropy-driven process

Multiscale Simulations of Biological Membranes : The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.Peer reviewe

Functionalization of the Parylene C Surface Enhances the Nucleation of Calcium Phosphate : Combined Experimental and Molecular Dynamics Simulations Approach

Interactions at the solid-body fluid interfaces play a vital role in bone tissue formation at the implant surface. In this study, fully atomistic molecular dynamics (MD) simulations were performed to investigate interactions between the physiological components of body fluids (Ca2+, HPO42-, H2PO4-, Na+, Cl-, and H2O) and functionalized parylene C surface. In comparison to the native parylene C (-Cl surface groups), the introduction of -OH, -CHO, and -COOH surface groups significantly enhances the interactions between body fluid ions and the polymeric surface. The experimentally observed formation of calcium phosphate nanocrystals is discussed in terms of MD simulations of the calcium phosphate clustering. Surface functional groups promote the clustering of calcium and phosphate ions in the following order: -OH > -CHO > -Cl (parent parylene C) approximate to -COO-. This promoting role of surface functional groups is explained as stimulating the number of Ca2+ and HPO42- surface contacts as well as ion chemisorption. The molecular mechanism of calcium phosphate cluster formation at the functionalized parylene C surface is proposed.Peer reviewe

Physiologically-relevant levels of sphingomyelin, but not GM1, induces a beta-sheet-rich structure in the amyloid-beta(1-42) monomer

To resolve the contribution of ceramide-containing lipids to the aggregation of the amyloid-β protein into β-sheet rich toxic oligomers, we employed molecular dynamics simulations to study the effect of cholesterol-containing bilayers comprised of POPC (70% POPC, and 30% cholesterol) and physiologically relevant concentrations of sphingomyelin (SM) (30% SM, 40% POPC, and 30% cholesterol), and the GM1 ganglioside (5% GM1, 70% POPC, and 25% cholesterol). The increased bilayer rigidity provided by SM (and to a lesser degree, GM1) reduced the interactions between the SM-enriched bilayer and the N-terminus of Aβ42 (and also residues Ser26, Asn27, and Lys28), which facilitated the formation of a β-sheet in the normally disordered N-terminal region. Aβ42 remained anchored to the SM-enriched bilayer through hydrogen bonds with the side chain of Arg5. With β-sheets in the at the N and C termini, the structure of Aβ42 in the sphingomyelin-enriched bilayer most resembles β-sheet-rich structures found in higher-ordered Aβ fibrils. Conversely, when bound to a bilayer comprised of 5% GM1, the conformation remained similar to that observed in the absence of GM1, with Aβ42 only making contact with one or two GM1 molecules. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.Peer reviewe

Dimerization of the pulmonary surfactant protein C in a membrane environment

Surfactant protein C (SP-C) has several functions in pulmonary surfactant. These include the transfer of lipids between different membrane structures, a role in surfactant recycling and homeostasis, and involvement in modulation of the innate defense system. Despite these important functions, the structures of functional SP-C complexes have remained unclear. SP-C is known to exist as a primarily alpha-helical structure with an apparently unstructured N-terminal region, yet there is recent evidence that the functions of SP-C could be associated with the formation of SP-C dimers and higher oligomers. In this work, we used molecular dynamics simulations, two-dimensional umbrella sampling, and well-tempered metadynamics to study the details of SP-C dimerization. The results suggest that SP-C dimerizes in pulmonary surfactant membranes, forming dimers of different topologies. The simulations identified a dimerization motif region V(21)xxxVxxxGxxxM(33) that is much larger than the putative A(30)xxxG(34) motif that is commonly assumed to control the dimerization of some alpha-helical transmembrane domains. The results provide a stronger basis for elucidating how SP-C functions in concert with other surfactant proteins.Peer reviewe

Maturation of the SARS-CoV-2 virus is regulated by dimerization of its main protease

SARS-CoV-2 main protease (Mpro) involved in COVID-19 is required for maturation of the virus and infection of host cells. The key question is how to block the activity of Mpro. By combining atomistic simulations with machine learning, we found that the enzyme regulates its own activity by a collective allosteric mechanism that involves dimerization and binding of a single substrate. At the core of the collective mechanism is the coupling between the catalytic site residues, H41 and C145, which direct the activity of Mpro dimer, and two salt bridges formed between R4 and E290 at the dimer interface. If these salt bridges are mutated, the activity of Mpro is blocked. The results suggest that dimerization of main proteases is a general mechanism to foster coronavirus proliferation, and propose a robust drug-based strategy that does not depend on the frequently mutating spike proteins at the viral envelope used to develop vaccines. (c) 2022 The Authors. Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).Peer reviewe

Reduced level of docosahexaenoic acid shifts GPCR neuroreceptors to less ordered membrane regions

G protein-coupled receptors (GPCRs) control cellular signaling and responses. Many of these GPCRs are modulated by cholesterol and polyunsaturated fatty acids (PUFAs) which have been shown to co-exist with saturated lipids in ordered membrane domains. However, the lipid compositions of such domains extracted from the brain cortex tissue of individuals suffering from GPCR-associated neurological disorders show drastically lowered levels of PUFAs. Here, using free energy techniques and multiscale simulations of numerous membrane proteins, we show that the presence of the PUFA DHA helps helical multi-pass proteins such as GPCRs partition into ordered membrane domains. The mechanism is based on hybrid lipids, whose PUFA chains coat the rough protein surface, while the saturated chains face the raft environment, thus minimizing perturbations therein. Our findings suggest that the reduction of GPCR partitioning to their native ordered environments due to PUFA depletion might affect the function of these receptors in numerous neurodegenerative diseases, where the membrane PUFA levels in the brain are decreased. We hope that this work inspires experimental studies on the connection between membrane PUFA levels and GPCR signaling. Author summary Our current picture of cellular membranes depicts them as laterally heterogeneous sheets of lipids crowded with membrane proteins. These proteins often require a specific lipid environment to efficiently perform their functions. Certain neuroreceptor proteins are regulated by membrane cholesterol that is considered to be enriched in ordered membrane domains. In the brain, these very same domains also contain a fair amount of polyunsaturated fatty acids (PUFAs) that have also been discovered to interact favorably with many receptor proteins. However, certain neurological diseasesassociated with the inadequate functioning of the neuroreceptorsseem to result in the decrease of brain PUFA levels. We hypothesized that this decrease in PUFA levels somehow inhibits receptor partitioning to cholesterol-rich domains, which could further compromise their function. We verified our hypothesis by an extensive set of computer simulations. They demonstrated that the PUFA-receptor interaction indeed leads to favorable partitioning of the receptors in the cholesterol-rich ordered domains. Moreover, the underlying mechanism based on the shielding of the rough protein surface by the PUFAs seems to be exclusive for multi-helical protein structures, of which neuroreceptors are a prime example.Peer reviewe

The role of hydrophobic matching on transmembrane helix packing in cells

Folding and packing of membrane proteins are highly influenced by the lipidic component of the membrane. Here, we explore how the hydrophobic mismatch (the difference between the hydrophobic span of a transmembrane protein region and the hydrophobic thickness of the lipid membrane around the protein) influences transmembrane helix packing in a cellular environment. Using a ToxRED assay in Escherichia coli and a Bimolecular Fluorescent Complementation approach in human-derived cells complemented by atomistic molecular dynamics simulations we analyzed the dimerization of Glycophorin A derived transmembrane segments. We concluded that, biological membranes can accommodate transmembrane homo-dimers with a wide range of hydrophobic lengths. Hydrophobic mismatch and its effects on dimerization are found to be considerably weaker than those previously observed in model membranes, or under in vitro conditions, indicating that biological membranes (particularly eukaryotic membranes) can adapt to structural deformations through compensatory mechanisms that emerge from their complex structure and composition to alleviate membrane stress. Results based on atomistic simulations support this view, as they revealed that Glycophorin A dimers remain stable, despite of poor hydrophobic match, using mechanisms based on dimer tilting or local membrane thickness perturbations. Furthermore, hetero-dimers with large length disparity between their monomers are also tolerated in cells, and the conclusions that one can draw are essentially similar to those found with homo-dimers. However, large differences between transmembrane helices length hinder the monomer/dimer equilibrium, confirming that, the hydrophobic mismatch has, nonetheless, biologically relevant effects on helix packing in vivo.Peer reviewe